Spin bath-mediated decoherence in superconductors 
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We consider a SQUID tunneling between 2 nearly degen- 
erate flux states. Decoherence caused by paramagnetic and 
nuclear spins in the low-T limit is shown to be much stronger 
than that from electronic excitations. The decoherence time 
is determined by the linewidth E of spin bath states, 
which can be reduced by a correct choice of ring geometry 
and isotopic purification. E can be measured in either field 
sweep or microwave absorption experiments, allowing both a 
test of the theory and design control. 



The fundamental importance of decoherence in Na- 
ture has been emphasized in fields ranging from quantum 
gravity and measurement theory to quantum computing 
and mesoscopic physics. Serious analysis of the mecha- 
nisms controlling decoherence began with work on dissi- 
pative tunneling and on "Macroscopic Quantum Coher- 
ence" (MQC) in SQUIDs To see MQC in SQUIDs 
requires T,f,(T)\A\ ^> 2-kTi, where is a "decoherence 
time", and A the tunneling matrix element. However, 
despite the successful observation of macroscopic tunnel- 
ing (MQT) in SQUIDs §, MQC has not been found. 
The unexpected low-T saturation of (T) in conductors 
, indicates we may not understand decoherence even 
in metallic systems (let alone in "qubits" or quantum 
computers!). 

Our thesis herein is that at low T in SQUIDs, t^(T) 
is controlled not by electronic or "oscillator bath" envi- 
ronments but by a "spin bath" Q of localized 
modes, including nuclear and paramagnetic spins as well 
as charge defects. Although spin bath environments have 
received considerable attention in nanomagnets [ ^0|pd| 
(and now in mesoscopic conductors |6|) their discussion 
for SQUID tunneling has been only sporadic ]j|JI^|. Here 
we show how they strongly suppress MQC (but usually 
with rather weak effects on MQT). Our results should 
allow experimentalists to both parametrize the spin bath 
effects, and quantitatively test the theory. 

1. Effective Hamiltonian. Ignoring electronic and elec- 
tromagnetic dissipation mechanisms the "bare" 
DC SQUID Hamiltonian for flux tunneling is 
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where C is the junction capacitance, (f> — 27r<I>/<i>o 
with $o the flux quantum, = —ihd/dQ, and <p x = 
27r$a;/$0! where $> x is the externally applied flux. Then 
Uq = <f>o/47r 2 L and g = 2ttLI c /^ , for a ring induc- 
tance L and junction critical current I c . In MQC or 
qubit designs <f> x <C 1 and tunnels between the two 
lowest wells, centred at <j> = <j>± = ±<f> m . One also 
assumes ksT < 7if2o/27r, where ilo is the Josephson 
plasma frequency of small oscillations in these wells. The 
system then truncates |l],Q,[7| to a two-level Hamilto- 
nian H = £ (t)r z + At x , where f is a Pauli vector, 
£o(t) ~ Uo4> m <p x (t) is a bias which can be varied in time, 
and current designs have |A/7i| as high as ~ 1 GHz. 

Consider now a set of TV two-level systems {<Jk} = 
{■Sfc,2fc}, representing paramagnetic ({.?&}) and nuclear 
({Ik}) spins jig] , at positions {r^}. Those in the SQUID 
couple to the conduction electron spin density s con( j(r) 
via 
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Jkdk ■ Scond(rfc) + 7fc<Tfc ■ B(r k ) 



(2) 



Spins in the substrate couple to the flux via the 2nd 
term. The {jk} = gk[J-k are spin magnetic moments, 
and the { Jk} represent electronic exchange for the {sfc}, 
and hyperfine coupling for the {Ik}- All {sk} are as- 
sumed paramagnetic Jl4| ], with Kondo energies ksTx ~ 
(J£ f ) 1 / 2 e- 1 / JJV (°) < k B T, where N(Q) is the Fermi- 
level density of states. The field B(r) = B x + B$(t) + 
^ fc 6fe(r), where B x is the external field, B$(r) comes 
from the SQUID supercurrent, and 
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is the dipolar field from 07. if 
conducting penetration depth 
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^From H — H + H sp we now derive a low-energy 
Hamiltonian, using the usual instanton method 0] 
wherein tunneling of <f> occurs in imaginary time r. A 
transition at r = gives a variation </>(r) = 4>mf(j), 
with /(t) ~ ± tan. _1 (e noT ). Consider now the local fields 
^k = 7fc-B(rfe) acting on the which during the in- 

stanton evolve like: 



Wfc(r) = uj^mk + wj^(r) , 



(4) 
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(cf. Fig. 1), where \fhk\ = 1 and /fc(r) evolves from if 
to — if during the instanton; \lf\ = 1 at the end-points. 
Then, defining B^(r k ) = [£?0 =+0m (r fc ) ± B 0= _ 0m (r fc )], 

^™fc = 7fc{Sx + J2 + \s+{v k )} (5) 

3 

4 = lk\\B^{v k )\ = l - lk 5B k (6) 

The time varying uAlk(r) causes transitions of a k ; writing 
\af >= Tjf\crZ >, which defines a transition matrix Tjf 
for the passage 17 — > if , we have in general 

Tf = exp j-i y + dTwj|i";(r) ■<?*} = e -^+**-** , 

(7) 

operating on <7fc, where ip k ,a k are typically complex 
In the long interval (over times ~ ft/ A) between instan- 
tons, sits in a static field ujf = uj k (l k = if); but dur- 
ing instantons, Tf operates. This gives immediately the 
diagonal and non-diagonal terms in an effective Hamilto- 
nian of "Central Spin" form §,0, 

valid at energy scales 

< hn -. 

H cff = [Af+e~ i £* dW ' i +H.c. 

+ ^2 {^k^k ■ rh k + f z u\a k ■ Ik) + io{t)r z , (8) 
k 

where the (p k are absorbed into the physical A. Eigen- 
states \a) of t z , with u = ±1 for t z =f , |, are converted to 
a pair of 2 iV -fold multiplets of coupled SQUID/spin bath 
states, with linewidth Eq = ^ k {w\) 2 and normalised dis- 
tributions w a (0 = (2/7r£;2)i/2 e -2(e- CT e ) 2 /£; 2 ; m energy 

£. The parameter E will be crucial to decoherence- note 
it only depends on the change SB k in field in (||), and is 
much easier to determine than the total field B(r k ). 

The above derivation ignores the slow spin diffusion 
and spin-lattice relaxation in the spin bath (which cause 
the vectors l k ,fn k to become dynamic variables). This 
occurs on timescales ~ (is or longer, ie., 3> ft/A, and is 
not relevant to decoherence in the present problem. 

2. Example: The parameters uif, u>\ are determined 
from the fields B^(r k ) (ie., knowing the supercurrent dis- 
tributions corresponding to <f> — ±4> m ), once we know 
{jk} and for all relevant nuclear and paramagnetic 
spins. We therefore assume homogeneous concentrations 
Xy j X J j and x s of paramagnetic spins in the ring, junction, 
and substrate respectively, for the geometry in Fig. 2, 
as well as a single nuclear species, with one nucleus per 
unit cell. Simple magnetostatics gives the results in Ta- 
ble I for ujf and cjjj. The vector a k depends on the 
detailed path followed by Zfc(r) during tunneling. Since 



Zfc(r) changes very quickly (on a timescale ^ D 1 ), from 
time-dependent perturbation theory \a k \ ~ oj\/M1 , and 
the mean number of bath spins flipping per transition is 
A ~ | J2 k \^k\ 2 ~ -E"o/ft 2 ^o (with an associated decoher- 
ence rate A A/ft ~ AE 2 /h 3 ttl). Thus A < 1 unless E a is 
unreasonably large. 

The dependence on sample size is the most striking 
result in Table I and Fig. 2, the contributions of all spins 
to E Q increasing rapidly as R decreases- the increase in 
the number of spins with R is more than offset by the 
weaker fields. The contribution from the substrate can 
easily be reduced by making it very thin and surrounding 
it with He, and one may also reduce E Q by reducing 
d,h,x, and cj> m Q. 

3. SQUID dynamics: A proper calculation of the 
SQUID dynamics requires adding to (||) a coupling to 
an oscillator bath, representing electrons, photons, and 
phonons. In the absence of a spin bath, the dom- 
inant electronic contribution gives a decoherence rate 
0,| r£ = IGir^kBT/Re 2 = & e kT/h where a s = 
{IQ-K^lfino/E^Q- 1 , Q(T) is the SQUID Q-f actor, and 
E c the junction charging energy. In a simple RCSJ model 
with shunt resistance R s and junction resistance Rj one 
has R- 1 = R^+RJ 1 (and Rj(T) ~ i? (I + e ABCS / T )/2, 
where Abcs is the superconducting gap and R a the 
normal junction resistance), so that Q and a s are T- 
independent at low T . 

Let us now calculate the SQUID dynamics for the 
Hamiltonian (ih, and thence the spin bath contribution 
to ; we then compare this to the electronic contribu- 
tion. Assuming A <C I (see above) we drop the {a k } from 
(||). Since u\, ojf <C A, we can treat these couplings per- 
turbatively; and a quick check of the numbers in Table I 
shows uj\ <C u)f almost always. Since in general rh and 
I are neither parallel nor perpendicular, we choose m as 
the spin quantization axis z. The component of I parallel 

to rfi is dealt with by redefining £ — > £+X)/c W I cos(/rn)cr^, 
which now depends on the spin bath state. Because l x 
has a transverse component, whenever the SQUID state 
changes the coupling term ^ fe ^\,S' k t z forces environ- 
mental spins to precess in a new local magnetic field, 
which can be viewed quantum mechanically as SQUID- 
induced transitions between the environmental states [|| . 
To quantify this effect we calculate the time correlation 
function P^(t), the probability JD that t z = 1 at time t 
(ie., <j) = +4> m ) if it was I at t = 0, after integrating out 
the spin bath. An instanton expansion gives (assuming 
e«|A|): 

2n 2m 

^Tt = |iA| 2( -" +m ' ) Y[ dt a Y[ dt' b F [r z (t), T z (t')] (9) 

nm a—1 b—1 

summed over "outgoing" and "return" paths T z (t),T z (t). 
The influence functional F p9| given from (@) , assuming 
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0, < fc B T, and uA < ^ 
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| / dseK»* s [T z (s)-T>(s)]\ 2 (10) 



(assuming an initial thermal environmental state). We 
distinguish 2 regimes: 

(a) Strong decoherence regime (E ^> A): Here F is 
negligible if instanton/anti-instanton pairs are separated 
by times > h/E Q , and we find immediately that F ~ 
exp[-£^(ti - t i+1 ) 2 /2?i 2 ], implying P TT (i) ~ (1/2)[1 + 
e-* /r «], with the rate 
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ie., completely incoherent quantum relaxation. 

(b) Weak decoherence regime (J5„ <C A): This problem 
is easily solved by going to the basis of eigenstates of f x . 
Their associated spin multiplets are widely separated in 
energy, so real transitions between them are impossible, 
and a perturbation expansion in w|/A gives 



H, 



off 



kk 1 



2A 



over 
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Tx + Yui*t (12) 



plus terms ~ 0({w\ ) 4 /A 3 ). Assuming £ < fc B T (so all 
states of a multiplet have equal thermal weight), we find 
that 



£44/8 A 
feft' 



E 2 J%hA 



(13) 



ie., roughly At^/2-kU ~ (A/E ) 2 oscillations survive be- 
fore phase decoherence sets in. We emphasize that no 
energy relaxation from |1) to |0) is involved here- this 
requires coupling to electronic excitations 

Comparing now electronic and spin bath decoherence 
rates, we see a crossover from electronic-dominated to 
spin bath-dominated decoherence around a temperature 
T c ~ (27r) 1 / 2 A 2 /fcs(5 s ii ; (strong decoherence regime), 
or T c ~ E 2 /8kBCt s A (weak decoherence regime). In the 
former case coherence will never be seen at any T; in the 
latter case the decoherence time will saturate at the value 
in (|l3|), below T c . In both cases the low T decoherence is 
controlled by the spin bath. We emphasize, on the other 
hand, that the spin bath will be almost invisible in MQT 
experiments- it causes litle dissipation, and adds to the 
tunneling exponent a factor 5S ~ nE o /%0. o , so typically 
SS <C 1 (and is T-independent!). 

4. Connection to Experiments: The concentration, 
type, and location of the paramagnetic impurities in the 
sample will usually be very uncertain. Luckily, ( |TT| ) and 
( |l2| ) show we only need to know Eq to parametrise the 
spin bath effects. In both regimes it can be determined by 



dynamic "fast passage" resonant tunneling experiments, 
from which one can extract W(e) = W+(e) + W-(e) by 
inverting the data pl)| . In the weak decoherence regime 
we can also look directly at the lineshape via microwave 
absorption between |0) and |1) manifolds. 

What this means is that (i) we may characterise the 
spin bath, extracting E Ql via well-established experimen- 
tal techniques, and then (ii) test the theory herein by 
comparing the predictions for P\-\{t) in (11) and ( fl3| ) 
with experiment. Note, incidentally, that just as in 
the nanomagnetic case p|-|Tl]], wide T-independent res- 
onant peaks in the sweep experiments are circumstantial 
evidence for a spin bath-mediated mechanism (oscilla- 
tor bath-mediated relaxation rates are T-dependent and 
typically increase as one moves further from resonance). 
Such peaks (of width ~ 0.4 K) were seen in recent ex- 
periments |Q, indicating a value E a ~ 0.2 K for this 
particular SQUID. If one is to see MQC, or to make su- 
perconducting qubits, E a must be reduced by at least 
10 2 . Ways to do this were indicated by our example- one 
wants very pure rings (including even isotopic purifica- 
tion if possible) with large R, small h, small junctions, 
thin substrates, and small <f> m |l6[| . 

In this paper we have shown how the low-T decoher- 
ence time in a SQUID must saturate at a value controlled 
by coupling to the spin bath. We thank the Institute 
for Theoretical Physics in Santa Barbara, where some 
of this work was done, and grant INTAS-2124 from the 
European community. 

Figures 

Fig. 1 (a) Evolution in imaginary time r of the field 
cuk acting on during tunneling of <f> from — <p m to 
+4> m - The field begins at uj^ and end at (2%. (b) A 
typical "path" (in the path integral sense) for r z (t), as 
the SQUID tunnels between | — <j) m > and \4> m >. The 
transition occurs in a "bounce time" 
between transitions ~ 2irn/ A. 



fl ; the time 



Fig. 2 A model DC SQUID ring, of height h, radius 
R, on a substrate of volume R 3 . The junction is a weak 
link of length d and radius r. We also show, for this 
example, the various contributions to E Q as a function of 
R, assuming h = 1 /im, r — 50 nm, d = 200 rim, and an 
insulating substrate with a concentration xs = 10~ 5 of 
paramagnetic impurities. We suppose ring and junction 
are made from Al (with j$ = 11.094M Hz/T, and I = 
5/2), with concentrations xj,x r — 10 -6 paramagnetic 
impurities. For definiteness we assume all paramgnetic 
impurities have spin s = 1, 7 S = 2fis = l&GHz/T '. 
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Ring Junction Substrate 

iVjv 4:Rh\Ln,N nr^dnN ~ BFun 

N s x.jN N xjN n x s N n 

SB k $ /i? 2 $o/(Rr) $o/R 2 

wf -fk8B k /2 7kSB k /2 -ykdBk/2 

Eg / 7 iv$o(AL/m J v) 1 / 2 /^ 3/2 J 7 jv$oMn JV ) 1 / 2 /2i? IlN^onJ, 2 /2R 1 ' 2 

E S Q s~f s $o(xj\ L hn N ) 1/2 /R 3/2 sj s ® (xjTTdn N ) 1/2 /2R s^ a $ (x a n N ) 1/2 /2R 1/2 

Un Jn[$o/R 2 + (^ /47r)(7jv + ^sXj,s)C N n N ] 

uj- 7 s [*o/-R 2 + (^o/47r)(7iv + ■ysXj l s)C N n N ] 

\dk\ "f k $o/2R 2 hn o ~{ k § /2RrhQ. lk $ /2R 2 hfl 

TABLE I. Parameters for the SQUID in Fig. 2, for nuclear (N) and paramagnetic (S) spins in the bulk ring, the junction, and 
the substrate. Nn and N s count all spins within a penetration depth Al of the surface - in the ring we assume h > \l > r (if 
Al > h, then substitute h instead of A). We assume B x ~ &o/R 2 - The number density njv of nuclear spins {I k } (with \I k \ = J) 
is ~ 1/fflo, where do is the lattice parameter; xj and x s are paramagnetic impurity concentrations per site (with \s k \ = s) in 
SQUID and substrate. Values for cj ± for ring and substrate spins (left blank in the Table) are the same as for the junction; 
Cat ~ 5 — 10 is a geometrical factor describing the effective number of nearest neighbour spins. Finally, E 2 = ^2 k (oJ k ) 2 , so 
E a ~ N} /2 luI and Eg ~ N^ 2 ^ (see text). 
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